Model of an exotic chiral superconducting phase in a graphene bilayer 
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We theoretically demonstrate the formation of a new type of unconventional superconductivity 
in graphene materials, which exhibits gapless property. The studied superconductivity is based 
on an interlayer pairing of chiral electrons in bilayer graphene, which results in an exotic s-wave 
spin-triplet condensate order with anomalous thermodynamic properties. These include the possi- 
bility of a temperature-induced condensation causing an increase of the pairing gap with increasing 
temperature, and an entropy of the stable superconducting state which can be higher than its value 
in the normal state. Our study reveals the analogy of the interlayer superconductivity in graphene 
materials to the color superconductivity in dense quark matter and the gapless pairing states in 
nuclear matter and ultra-cold atomic gases. 
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The microscopic Bardeen-Cooper-Schricffcr (BCS) 
theory of superconductivity is based on the forma- 
tion of the so called Cooper pairs between electrons with 
opposite spin and equal but opposite momentum, when a 
weak attractive interaction exists between the electrons. 
Despite first being applied specifically to explain metallic 
superconductivity, the BCS pairing was discovered later 
on to be a general effect which could appear in other ar- 
eas of physics (see Ref. for a recent overview). While 
in conventional metallic superconductors an energy gap 
is opened in the electronic excitation spectrum on the 
Fermi surface^, unconventional pairing states can be 
realized in a wide range of fermionic systems in which 
superconductivity coexists with the gapless spectrum of 
the normal (N) phase. A great interest has been devoted 
to the BCS pairing in composite systems of two (or more) 
types of fcrmions with different Fermi surfaces, in which 
despite a nonvanishing pairing gap (PG) or order param- 
eter, there appears no energy gap on the Fermi surface 
due to a separation of N and superconducting (S) phases 
in the momentum space [ij . Gapless color superconduc- 
tivity in dense quark matter at low temperatures [H| and 
its analogous gapless pairing states in nuclear matter Q 
and ultra-cold atomic gases with two types of atoms or 
identical atoms with different spin states in an external 
magnetic (see e. g. Q), are famous examples of this 
kind. 

The color superconductivity is the state of matter at 
highest densities and low temperatures which can be 
reached at the interior of dens neutron stars [Bj]. In this 
Letter we report on the theoretical study of interlayer 
electronic pairing in pure bilayer graphene, which ex- 
plores the potential of this carbon-based material to lend 
itself to an experimentally feasible study of the relativis- 
tic gapless superconductivity. This may also simulate 
the regime of color superconductivity. Graphene @, the 
two dimensional solid of carbon atoms with honeycomb 
lattice structure, has already been proven to exhibit a va- 
riety of compelling pseudo-relativistic and quantum elec- 
trodynamics phenomena such as Klein tunneling [i| , Zit- 




FIG. I: (a) Sketch of the bilayer graphene lattice structure 
with the interlayer chiral superconducting correlations (illus- 
trated by wavy lines) between electrons of different types of 
sublattices. Dirac cone bands of two monolayers of graphene 
on the valleys K\ and K' 2 (right panel), with two states of 
wave vector \k) and | — k) which are coherent superpositions 
of the sublattice pseudospin states |^4i,2) and \B\^)- The in- 
terlayer pairing between | ± k) is asymmetrically partial, since 
only \A\) and \B 2 ) components are paired, (b) The mean 
field phase diagram of the chiral superconductivity showing 
dependence on the interlayer hopping energy t± and the tem- 
perature T for a pure bilayer graphene (p = 0). The line of 
transition from N to S phases is shown in which dashed and 
solid parts indicate the first and the second order transitions, 
respectively. Color plot of the pairing gap (PG) is also pre- 
sented, (c) The resulted gapless S band structure along k y 
for undoped bilayer graphene and for t± — 0A3t. 



terbewugung and universal light absorption[ll|, to 
name a few. This is due to its unique zero-gap elec- 
tronic band structure with conically shaped conduction 
and valence bands touching each other at the corners of 
hexagonal first Brillouin zone, known as Dirac points (sec 
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Fig. UK)- An important aspect in graphene is the connec- 
tion between this specific band structure and the sublat- 
tice pseudospin degree of freedom, which makes electrons 
behave like two dimensional massless Dirac fermions with 
a pseudo-relativistic chiral property [12j. Here we find 



that in bilayer graphene, the interplay between the inter- 
layer pairing of electrons with the same sublattice chiral- 
ity and the asymmetric arrangement of the sublattices of 
the two layers results in a gapless superconducting which 
posses an exotic s-wave spin triplet symmetry of the order 
parameter and an anomalous N-S phase diagram, shown 
in Fig. [TJ). We further find that the entropy in the chiral 
S state can be higher than the entropy in N state. This 
finding is in striking contrast to the common behavior of 
the entropy upon a phase transition. 

We start with a qualitative description of the interlayer 
pairing in bilayer graphene, which reveals the key fea- 
tures of the chiral superconductivity and its anomalous 
properties. As is shown in Fig. QJi, in bilayer graphene 
two monolayers with A\ and B\ triangular sublattices in 
top layer and A2 and B2 triangular sublattices in bot- 
tom layer, are arranged according to Bernal stacking in 
which every A\ site of top layer lies directly above a B2 
site of bottom layer. In the absence of hopping between 
the layers(ii_ = 0) the electronic states of the layers are 
not mixed in N state. In the layer 1(2) the electronic 
states lie in two Dirac-cone bands on the two inequiva- 
lent -Ki(2) and K'\{2) valleys (Dirac points). In each layer, 
pseudospin chirality links the relative amplitude of elec- 
tron wave function (to be in the A or B sublattices) to 
the direction of its momentum. Each momentum state 
in Dirac cones of layers 1,2 |fci,2), is a coherent superpo- 
sition of the two sublattice pseudospin states lAi^) and 
\Bi,i) as 1^1,2) = ^1,2) + exp(iip kl 2 )\Bi^ 2 ), with the an- 
gles <pki,2 determining the directions of the wave vectors 
fci.2- The interlayer superconductivity will pair electrons 
in Al and B2 sublattices, but not those in A2 and Bl 
sublattices, which implies that the pairing between the 
states and |A^) is asymmetrically partial with re- 
spect to the pseudospin degree of freedom. Taking into 
account the chirality, together with the requirement that 
the paired electrons have equal but opposite momentum 
on the valleys K\ and K 2 , we obtain a rough qualitative 
picture of N and S states configuration in the momentum 
space. In each valley, only electrons of one sublattice be- 
comes S and the other sublattice remains in N state. This 
simultaneous existence of N and S phases can give rise 
to the gapless property of the chiral superconductivity 
for an undoped bilayer. Having a partial pairing at zero 
temperature T = 0, a thermal excitation can redistribute 
electrons in N component of a wave vector states in new 
momentum states of coherent superposition of A and B 
states, where they have potential for a partial pairing 
again. This corresponds to a temperature-induced con- 
densation process in which two primarily normal elec- 
trons are thermally excited to the coherent states, where 



they can now be partially paired and absorbed into the 
chiral superconducting condensate. By this process the 
pair amplitude can increase with increasing temperature. 

The above described qualitative picture can be sup- 
ported by a mean field calculation based on tight-binding 
theory of 7r electrons in graphene. In the absence of su- 
perconductivity, we consider the following Hamiltonian 
1 1 31 ] for bilayer graphene 

#0 = -H n t ,i, a - t Y (°<ti,<A,j> + h - c -) W 



where at,%,a (&£,«, CT ) and a\ icr (b\ i ) are the annihilation 
and creation operators of an electron in the ith unit cell 
in sublattice A(B) and layer £(= 1,2); a = ± denotes 
the spin state of the electron; n^i t „ is the corresponding 
on-site particle density operator. The intralayer nearest 
neighbor hopping energy t « 3eV determines the Fermi 
velocity in graphene as vp rj 10 6 m/s, and t± ~ OAeV 
(t±/t ~ 0.13) is the dominant interlayer hopping energy 
between the nearest neighbors Al and B2; the chemical 
potential /x can be controlled by gate voltages. The in- 
terlayer superconductivity can be introduced by adding 
the following potential to Hq 



Vx = -g. 



2,i,ff' f '2,i,er' 



(2) 



where g± is the S coupling energy [3 EH ■ Below, we will 
discuss the origin of this attractive potential. The inter- 
layer Coulomb repulsion, competing with this attractive 
potential, is expected to be rather weak because of the 
interlayer spacing between the electrons. We note that 
the potential ([2]) presents an on-site local interaction in 
the 2D plane of the bilayer, which requires the pairing 
to have an isotropic s-wave orbital symmetry. On the 
other hand the pairing is anti-symmetric with respect to 
the pseudospin (i.e. sublattice) degree of freedom. Pauli 
exclusion principle requires that the total wave function, 
composed of the product of the orbital, spin and pseu- 
dospin components, should be anti-symmetric under the 
exchange of electrons. Therefore we conclude that the in- 
teraction ([2]) allows only an spin-triplet pairing. We find 
that the interacting Hamiltonian H = Hq + V± can be 
decoupled by the following spin-triplet order parameter 



(3) 



which describes the interlayer pairing of electrons with 
the sublattices Ai-B 2 , the valleys Ki-K 2 (and K' 1 -K 2 ), 
and the spins t~4- (and . It is worth to note that while 
the unconventional s-wave spin-triplet pairing is expected 
to exhibit an odd- frequency superconductivity [16(, we 
have found that the order parameter §3§ has an even 
parity in the frequency space. This is the result of an 
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FIG. 2: Exotic thermodynamic properties of the chiral su- 
perconductivity, (a) Pairing gap versus T for tj_/t = 
0.0, 0.13, 0.2, 0.34 (from top to bottom, respectively), (c) Dif- 
ference in entropies of N and S states versus T for the same 
values of parameters as in (a). 



FIG. 3: Stability of the exotic chiral superconducting states. 
Difference in thermodynamic potentials of N and S states ver- 
sus the pairing gap calculated at different temperatures for 
interlayer hopping energies O.lt (a) and 0.18t (b). 



odd-pscudospin parity, as described above. We obtain 
the excitation spectrum by replacing the mean field or- 
der parameter ((3]) in H = H + V±_ and performing the 
Fourier transformation of the Hamiltonian in the space 
of 2D wave vector k of electrons. For zero doping \i = 0, 
we obtain a simple expression for the spectrum E$} which 
reads 



a 



[^(Ax + 7t2j* + 4e2 + l(A ± + st±)], (4) 



where I, s,a = ± indicate different branches of the spec- 
trum. Here ek = ±f|7k| with 7^ = ^ 5 exp — ik.<5, 
with S indicating the nearest neighbor position vectors of 
graphene. The spectrum (j4]) represents a gapless super- 
conductivity as can be seen in Fig. [TJ;, in which different 
branches of Ef~ versus k y are plotted when fi = 0. Now, 
having the excitation spectrum and using the standard 
many-body Green's function method, we can obtain the 
self-consistent equation for the pairing gap (PG) Aj_, as 
follows 
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We can also obtain the thermodynamic potential fig and 
the entropy Ss of S state from the following formulas, 
respectively, 



n s = 



Ai 
9x 
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S S = -2k B In/£' + (1 - /^)ln(l - /£')], (7) 

k,s,l 

where f± = [1 + exp^-E^ is the Fermi-Dirac oc- 
cupation factor with j3 = 1/ksT (ks being Boltzmann 
constant). The expressions for the thermodynamic po- 
tential and entropy in N state follow from Eqs. © and 
by taking A_l = 0. 



The temperature dependence of PG for different values 
of the interlayer hopping energy can be obtained from a 
numerical solution of Eq. (|S|). We have examined the 
stability of the resulted S states by ensuring that any 
solution of §5§ presents indeed a global minimum of fis 
(see Fig. [3]). The resulting phase diagram and Aj_(T) 
is surprisingly anomalous, as is presented in Fig. [T}d and 
Fig. [2^, respectively, for an undoped bilayer (fj, = 0). 
For t± = 0, the gap increases from its zero temperature 
value with increasing temperature and reaches a maxi- 
mum before showing a BCS like decrease with T. As 
described above, this is the result of thermal excitations 
of electrons from a N valence band into a S conduction 
band, where they can form a condensate and contribute 
to the PG. An interlayer hoping t± introduces a nor- 
mal coupling between the Al and B2 sublattices and will 
weaken the superconducting correlation and hence the 
interlayer PG. The interplay between this effect and the 
temperature-induced increase of the chiral condensation 
leads to the most amazing behavior of PG. As is seen 
in Fig. [Tb and Fig. [2k., it appears that for t± exceeding 
some critical value, PG vanishes at a temperature range 
starting from T = 0, while it takes a sizable value at 
higher temperatures. 

In the phase diagram of the chiral superconductivity 
in t± — T plane (Fig. QJi) the color plot of PG is also 
shown. The normal-superconducting (NS) transition is 
first order at lower temperatures (the dashed part of a 
NS boundary), and second order at higher temperatures 
(the solid part). This can be seen from Fig. [3j in which 
the potential difference fts — Ojv in N and S states is 
plotted versus A^ for ^ = and two values of t±/t = 0.1 
and 0.18. While for t±/t = 0.1, the low temperature 
NS transition is continuous, for t±/t = 0.18, it happens 
discontinuously, as the stable solution Aj^ undergoes a 
jump from zero to a finite value at the transition. 

The anomalous behavior of A±(T) is associated with 
a very unusual behavior of the entropy, as is presented 
in Fig. by plotting the entropy difference Ss — Sn in 
N and S states as a function of T for different values of 
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t± . The entropy difference undergoes a change of sign at 
the temperature of maximum PG, which implies that the 
entropy of the S state is higher than that of the N state, 
in the temperature interval where PG increases with T. 
This unusual finding can also be understood in terms of 
the above described temperature-induced condensation 
process. Similar to the above found exotic behavior of 
PG and the entropy are also predicted to exists in the 
pairing state of asymmetric nuclear matter j(| , and gap- 
less color superconducting phase of dens quark matter 
in which correlations between quarks with three different 
colors and three different flavors results in a variety of 
pairing patterns [H, [I?} . A similar to colored pairing of 
quarks can be simulated by the interlayer chiral super- 
conductivity in trilayer graphene with two different ABA 
and ABC stackings [18J, which provides pairing between 
electrons whose states lie in three distinct Dirac cones 
with opposite sublattice pseudospins. 

Finally, let us comment on the experimental realiza- 
tion of our predicted exotic superconductivity. The pos- 
sibility of an intrinsic superconductivity, with plasmon 
or phonon mediated pairing interactions, has been pro- 



posed in graphene coated with a metal [15| . For a coated 
bilayer graphene the formation of the S state is expected 
to be closely similar to that of the graphite intercalated 
with alkaline metals, for which a critical temperature up 
to 11.5 K has been reported (see e. g. [13] )• This can 
be served as an estimation of the critical temperature 
(and essentially the coupling energy g±) for the inter- 
layer coupling. Alternatively, the superconductivity can 
be induced by the proximity to a metallic superconduc- 
tor. Already, significant progress has been made in fab- 
ricating highly transparent contacts between a graphene 
monolayer and a superconductor, TiAl bilayer is proved 
to be a suitable material because of good electrical con- 
tact of Ti to graphene (see e. g. [lj| ) . For the proximity 
induced superconductivity T c is determined by that of the 
contacted superconductor (for TiAl bilayer it is about 1.3 
K). A superconductor-bilayer-graphene-superconductor 
setup 2l|, with several S electrodes contacting indepen- 
dently to the top and bottom layers can provide the possi- 
bility for an induction of controllable interlayer supercon- 
ducting correlations. We anticipate that a hybrid system 
(like a SNS Josephson junction) made of chiral super- 
conductors to exhibit peculiar properties due to the ex- 
otic order parameter and thermodynamic characteristics 
of these type of superconductors. This can be used for 
the experimental observation of chiral superconductivity, 
which may also provide potential for new applications. It 
is worth to note that in both of the above proposed meth- 
ods, measurements sensitive to interlayer pairing, similar 
to those of the intrinsic Josephson effect in layered high- 
T c superconductors (see e. g. [23]) , are required. 

In conclusion, we have predicted a new type of gapless 
superconductivity which can be realized by an interlayer 
pairing of rclativistic-likc chiral electrons in graphene ma- 



terials. This chiral superconductivity posses an exotic s- 
wave spin-triplet order parameter with anomalous ther- 
modynamic properties. We have found that the pairing 
gap can increase with increasing of temperature, with 
the possibility of vanishing at a lower temperature range 
including T = 0, while having a finite value at higher 
temperatures. The entropy of these anomalous stable 
superconducting states is found to be higher than that of 
the normal state. 
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